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RTCC REQUIREMENTS FOR MISSIONS F AND G: 
SELENOGRAFHIC -MNBY VELOCITY TRANSFORMATION FORMULATION 

By Paul F. Flanagan 


SUMMARY 


This document presents the formulation and implementation procedures 
for use in the RTCC for transforming velocity vectors from the seleno- 
graphic coordinate system to the mean nearest Besselian year (MNBY) 
system. The formulation includes equations for computing the lunar 
angular velocity vector. 


INTRODUCTION 


The selenocentric MNBY system is moon centered and has its axes 
parallel to the following: the X-axis is in the mean earth equatorial 

plane and through the mean vernal equinox at the beginning of the nearest 
Besselian year; the Z-axis is north through the mean earth spin axis; 
and the Y-axis is on the mean earth equatorial plane and forms a right- 
handed orthogonal system. 

The moon fixed selenographic system is defined as follows. The 
X-axis is fixed in the moon's equatorial plane and points approximately 
toward the center of the earth; the Z-axis is fixed in the moon through 
the moon spin axis; and the Y-axis is in the equatorial plane and com- 
pletes the right-handed orthogonal system. 

The interpolated precession-libration matrix (L) used to transform 
the velocity is also used to transform position vectors. If 

position vector in the selenocentric MNBY system and R is a 
vector in the selenographic system, then 

r m = L \ 


*M ls 8 

position 


and 


R G = IJi M 
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VELOCITY VECTOR TRANSFORMATION 


The inertial MNBY velocity, R,., of a vehicle with a fixed sele- 
nographic position, R^ , is computed using the following equation 

u 

^ * *\j> 

where W is the moon's angular velocity expressed in the selenographic 
* system. 

The following general equations can be used to compute the inertial 
velocity for a vehicle moving relative to the selenographic system at 

velocity R„ or for the reverse computation, 
o 

\ • lT ^g + w x V 

*o c ^ - w * R 0 

COMPUTATION OF W 


The following equations are used to compute the angular velocity 
vector, W. These equations account for the rate of change of £1, the 
longitude of the mean ascending node of the lunar orbit measured in the 
ecliptic plane from the mean equinox of date, and M, the geocentric mean 
longitude of the moon measured in the ecliptic plane from the mean equinox 
of date to the mean ascending node of the lunar orbit, and then along the 
orbit. 

The computations neglect the small effect of the rates of precession, 
obliquity of the ecliptic, and physical libration. The inclination (i) 
of the lunar orbit to the ecliptic is constant and is 1 0 32'39" (ref. 1). 

Prelaunch initialization: 

1. Input year of launch, YOL, and day of launch, DOL. 

Ti .r.nary 1 is day 1. 

2. Compute Julian centuries from January 0.5, 1900 to midnight 
prior to launch: 


TIL = integral part [(YOL - 190l)/L] 

TB = (YOL - 1900)365 + TIL + DOL - 0.5 
TM = TB/36525 

• • 

3. Compute time rates of angles, Q and M in rad/hr: 
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n « [- 1931 *. 112008 3 + 0 . 00 U 155556 TM + ( 0.6666667 x io’ 5 )tm ]rate 

M *= LU81267.8831U17 - 0.002266667TM + (O .5666667 X 10 > )TI^]RATE 

where RATE = */(l80 • 36525 • 24) « 1.9110212776572292 x 10* and converts 

• • 

Q and M from deg/Julian century to rad /hr. 

1. Compute constant terms: 

K1 = f) cos i 
K2 = n sin i 

5* Compute in raa/hr: 

w z = M - n + K1 

Computation of and at any time t during mission: 

1. Input hours from midnight prior launch, H. 

2. Compute Julian centuries from January 0.5, 1900 to anytime tj 

T = (TB + H/21i)/ 36525 
5. Compute M - 0 in deg: 

M - 0 = 11.2508889 + 1*8 3202. 02515T 

- 0.003211111T 2 - (0.0333333 X 10*^)T 3 

1*. Compute W and W in rad/hr: 

X Y 

W = K2 sin (M - Q) 

W = K2 cos (M - n) 

A listing of a single-precision checkout program to compute the 
angular velocity, W, is given in table I. A sample output for RTCC 
checkout is provided in table II for midnight January 1 through 10, 1969* 


REPRODUCIBILITY OF THE ORIGINAL PAGE IS POOR. 


TABLE I.- LISTING OF SINGLE-PRECISION CHECKOUT PROGRAM 
TO COMPUTE THE ANGULAR VELOCITY VECTOR W 


r, "-1. r 

rT -f-Tf .(3) 

A r ci ; T[voi.,h 

i Apr. PTrrou 
Tvr’.rvoi. 

Tll.r ( Tvn| - j *>( l > //. 

T .r (YfL-r; . f . ) * y > f , -**T I ’ + r -11 

7 , z T ' / 3 7 s - ^ A . 

r r r .= 7 ),i.ur^L^/i c. 

Arfil.I = r *5.';P.*r ^A/5 >co. 

rrz^^.ur.cr^c.oo^ss^ »+.cs< < ?‘;«7 7 ) 7 ,4 )*r\ * * 

: ~ = (AM?77.sp31417M-.r0f;25S^+.f ) • T )-r- ' 

i?i = rr*r ■:[ tr cLi j 

!" r nr *f T * [ A f‘CL J 1 
•♦.C3) = (; r-nr-MM) 

Tz ( T?>+F/£4 •)/3^52f>. 

• v 0= C 1 1 .250RS°+( 4F3202. 025 1 5+(-. 00321 1 1 1 1 - . 0333333E-5*T ) *T ) *T ) *r>E oc a 
’J( l >z(K?*SINfrM0)) 

J(2)t(K?*con[fiMO]) 

77=1 ,0/( 36525,0*24.0) 

"( I )=W(| )*77 
"r?.)r ,t (P)*77 
"( 3) r <( 3) *77 
7 J 7H 1 ) f ’v(7) t 1(3) ) 

niqpLAV[voi.,roi. f ! i 

7 r T c i 

rr^ 


t 
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TABLE II.- ANGULAR VELOCITY VECTORS FOR 
PRELAUNCH INITIALIZATION FOR MIDNIGHT OF JANUARY 1 THROUGH 10, 1969 


Angular velocity, rad/hr 

W 

w 

w 

X 

y 

z 

-o. 9291^761 U x I0" fj 

-0.U5U106239 x io“ G 

0.958135087 X 10“ 2 

-0.1008141172 x icf 5 

-0.229I4200U1 x io“ G 

0.958135087 x 10“ 2 

-0.103U15296 X lO" 5 

0.7^4292239 x io' H 

0.958135087 X 10“ 2 

-0.100500508 X 10“ 6 

0.2U 391089 x 10“ 6 

0.958135087 x 10’ 2 

-0.922515165 * io“ 6 

0.467432877 x 10“ G 

0.958135087 X 10“ 2 

-0.791061U5I x 10“ G 

0.666145255 X 10“ 6 

0.958135037 X 10“ 2 

-0.617621078 x 10“ & 

0.829501023 X ]0” G 

0.958135087 x 10“ 2 

-0.U11399568 X 10“ 6 

O.9I48829879 X 10” G 

0.958135087 X 10" 2 

-0.1833U2UU8 X i0“ e 

0.101779825 X 10“ r> 

0.958135087 X 10“ 2 

0.5^^^580147 X 10” 7 

O.IO327I4556 X 10~ b 

0.958135087 X 10“ 2 
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